32 research outputs found
Recommended from our members
New Direction in Hydrogeochemical Transport Modeling: Incorporating Multiple Kinetic and Equilibrium Reaction Pathways
At least two distinct kinds of hydrogeochemical models have evolved historically for use in analyzing contaminant transport, but each has important limitations. One kind, focusing on organic contaminants, treats biodegradation reactions as parts of relatively simple kinetic reaction networks with no or limited coupling to aqueous and surface complexation and mineral dissolution/precipitation reactions. A second kind, evolving out of the speciation and reaction path codes, is capable of handling a comprehensive suite of multicomponent complexation (aqueous and surface) and mineral precipitation and dissolution reactions, but has not been able to treat reaction networks characterized by partial redox disequilibrium and multiple kinetic pathways. More recently, various investigators have begun to consider biodegradation reactions in the context of comprehensive equilibrium and kinetic reaction networks (e.g. Hunter et al. 1998, Mayer 1999). Here we explore two examples of multiple equilibrium and kinetic reaction pathways using the reactive transport code GIMRT98 (Steefel, in prep.): (1) a computational example involving the generation of acid mine drainage due to oxidation of pyrite, and (2) a computational/field example where the rates of chlorinated VOC degradation are linked to the rates of major redox processes occurring in organic-rich wetland sediments overlying a contaminated aerobic aquifer
Recommended from our members
Wavelet-based local mesh refinement for rainfall–runoff simulations
A wavelet-based local mesh refinement (wLMR) strategy is designed to generate multiresolution and unstructured triangular meshes from real digital elevation model (DEM) data for efficient hydrological simulations at the catchment scale. The wLMR strategy is studied considering slope- and curvature-based refinement criteria to analyze DEM inputs: the slope-based criterion uses bed elevation data as input to the wLMR strategy, whereas the curvature-based criterion feeds the bed slope data into it. The performance of the wLMR meshes generated by these two criteria is compared for hydrological simulations; first, using three analytical tests with the systematic variation in topography types and then by reproducing laboratory- and real-scale case studies. The bed elevation on the wLMR meshes and their simulation results are compared relative to those achieved on the finest uniform mesh. Analytical tests show that the slope- and curvature-based criteria are equally effective with the wLMR strategy, and that it is easier to decide which criterion to take in relation to the (regular) shape of the topography. For the realistic case studies: (i) slope analysis provides a better metric to assess the correlation of a wLMR mesh to the fine uniform mesh and (ii) both criteria predict outlet hydrographs with a close predictive accuracy to that on the uniform mesh, but the curvature-based criterion is found to slightly better capture the channeling patterns of real DEM data
Evaluation of the field-scale cation exchange capacity of Hanford sediments
Three-dimensional simulations of unsaturated flow, transport, and multi-component, multi-site cation exchange in the vadose zone were used to analyze the migration of a plume resulting from a leak of the SX-115 tank at the Hanford site, USA. The match within about 0.5 meters of the positions of retarded sodium and potassium fronts suggests that the laboratory-derived parameters may be used in field-scale simulations of radionuclide migration at the Hanford site
Recommended from our members
Kaolinite dissolution and precipitation kinetics at 22oC and pH 4
Dissolution and precipitation rates of low defect Georgia kaolinite (KGa-1b) as a function of Gibbs free energy of reaction (or reaction affinity) were measured at 22 C and pH 4 in continuously stirred flowthrough reactors. Steady state dissolution experiments showed slightly incongruent dissolution, with a Si/Al ratio of about 1.12 that is attributed to the re-adsorption of Al on to the kaolinite surface. No inhibition of the kaolinite dissolution rate was apparent when dissolved aluminum was varied from 0 and 60 {micro}M. The relationship between dissolution rates and the reaction affinity can be described well by a Transition State Theory (TST) rate formulation with a Temkin coefficient of 2 R{sub diss} (mol/m{sup 2}s) = 1.15 x 10{sup -13} [1-exp(-{Delta}G/2RT)]. Stopping of flow in a close to equilibrium dissolution experiment yielded a solubility constant for kaolinite at 22 C of 10{sup 7.57}. Experiments on the precipitation kinetics of kaolinite showed a more complex behavior. One conducted using kaolinite seed that had previously undergone extensive dissolution under far from equilibrium conditions for 5 months showed a quasi-steady state precipitation rate for 105 hours that was compatible with the TST expression above. After this initial period, however, precipitation rates decreased by an order of magnitude, and like other precipitation experiments conducted at higher supersaturation and without kaolinite seed subjected to extensive prior dissolution, could not be described with the TST law. The initial quasi-steady state rate is interpreted as growth on activated sites created by the dissolution process, but this reversible growth mechanism could not be maintained once these sites were filled. Long-term precipitation rates showed a linear dependence on solution saturation state that is generally consistent with a two dimensional nucleation growth mechanism following the equation R{sub ppt}(mol/m{sup 2}s) = 3.38 x 10{sup -14} exp[- 181776/T{sup 2} 1n{Omega}]. Further analysis using Synchrotron Scanning Transmission X-ray Microscopy (STXM) in Total Electron Yield (TEY) mode of the material from the precipitation experiments showed spectra for newly precipitated material compatible with kaolinite. An idealized set of reactive transport simulations of the chemical weathering of albite to kaolinite using rate laws from HELLMANN and TISSERAND (2006) and this study respectively indicate that while pore waters are likely to be close to equilibrium with respect to kaolinite at pH 4, significant kaolinite supersaturation may occur at higher pH if its precipitation rate is pH dependent
Recommended from our members
Mineral dissolution kinetics at the pore scale
Mineral dissolution rates in the field have been reported to be orders of magnitude slower than those measured in the laboratory, an unresolved discrepancy that severely limits our ability to develop scientifically defensible predictive or even interpretive models for many geochemical processes in the earth and environmental sciences. One suggestion links this discrepancy to the role of physical and chemical heterogeneities typically found in subsurface soils and aquifers in producing scale-dependent rates where concentration gradients develop. In this paper, we examine the possibility that scale-dependent mineral dissolution rates can develop even at the single pore and fracture scale, the smallest and most fundamental building block of porous media. To do so, we develop two models to analyze mineral dissolution kinetics at the single pore scale: (1) a Poiseuille Flow model that applies laboratory-measured dissolution kinetics at the pore or fracture wall and couples this to a rigorous treatment of both advective and diffusive transport, and (2) a Well-Mixed Reactor model that assumes complete mixing within the pore, while maintaining the same reactive surface area, average flow rate, and geometry as the Poiseuille Flow model. For a fracture, a 1D Plug Flow Reactor model is considered in addition to quantify the effects of longitudinal versus transverse mixing. The comparison of averaged dissolution rates under various conditions of flow, pore size, and fracture length from the three models is used as a means to quantify the extent to which concentration gradients at the single pore and fracture scale can develop and render rates scale-dependent. Three important minerals that dissolve at widely different rates, calcite, plagioclase, and iron hydroxide, are considered. The modeling indicates that rate discrepancies arise primarily where concentration gradients develop due to comparable rates of reaction and advective transport, and incomplete mixing via molecular diffusion. The magnitude of the reaction rate is important, since it is found that scaling effects (and thus rate discrepancies) are negligible at the single pore and fracture scale for plagioclase and iron hydroxide because of the slow rate at which they dissolve. In the case of calcite, where dissolution rates are rapid, scaling effects can develop at high flow rates from 0.1 cm/s to 1000 cm/s and for fracture lengths less than 1 cm. At more normal flow rates, however, mixing via molecular diffusion is effective in homogenizing the concentration field, thus eliminating any discrepancies between the Poiseuille Flow and the Well-Mixed Reactor model. This suggests that a scale dependence to mineral dissolution rates is unlikely at the single pore or fracture scale under normal geological/hydrologic conditions, implying that the discrepancy between laboratory and field rates must be attributed to other factors
Long-term fluxes of reactive species in macrotidal estuaries: estimates from a fully transient, multicomponent reaction-transport model
A coupled, fully transient, multicomponent reaction-transport model has been developed to estimate long-term fluxes of reactive compounds in strong tidal estuaries. The model is applied to a preliminary analysis of the carbon cycle in the Scheldt estuary in Belgium and The Netherlands. The model provides a realistic description of the residual circulation in a strong tidal estuary and includes the essential feedback mechanisms between interdependent chemical species. The model has been used to analyze the fundamentally transient nature of strong tidal estuaries and, in particular, the effect of these non-steady state conditions on the long-term fluxes of chemical species out of the estuary. The results indicate that flux estimation techniques based upon steady-state assumptions may result in significant errors. The model has also been used to investigate biogeochemical interactions characterized by a large spectrum of time scales, which it does by including simultaneous equilibrium reactions and kinetically-mediated processes. Simulations carried out with the model suggest that a formulation based upon microbially-mediated, kinetically-controlled reactions provides a superior description of solute profiles in the Scheldt estuary than does a global equilibrium redox formulation. The mixed equilibrium-kinetic formulation also makes it possible to track simultaneously two master variables: the redox state of the system and the pH. By providing strong constraints on the system, these two master variables can be used to test the model's self-consistency. The simulations carried out with the model suggest the pH profile in the Scheldt estuary is the result of a balance of biogeochemical reactions which produce H+ and degassing which consumes H+ and not the result of simple mixing between seawater and freshwater